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We consider a one dimensional spin 1/2 chain with Heisenberg interaction in a disordered parallel 
magnetic field. This system is known to exhibit the many body localization (MBL) transition at 
critical strength of disorder. We analyze the response of the chain when additional perpendicular 
magnetic field is applied to an individual spin and propose a method for accurate determination 
of the mobility edge via local spin measurements. We further demonstrate that the exponential 
decrease of the spin response with the distance between perturbed spin and measured spin can be 
used to characterize the localization length in the MBL phase. 


Introduction.- Theoretical studies of localization in 
many body systems were originally focused on interacting 
electrons in disordered metals [1-3], where experimental 
observation of localization remains a challenging task due 
to electron-phonon interaction that ruins many-particle 
states by introducing external decoherence to the system 
of interacting electrons in real metals. The interest in 
observation of many-body localization (MBL) has shifted 
to artificial quantum systems well isolated from their en¬ 
vironment such as ultra-cold atomic gases [4], trapped 
ions [5, 6] and superconducting circuits containing many 
interacting qubits [7]. These systems have configurable 
Hamiltonians and their observables can be measured with 
high precision, which gives another advantage to artifi¬ 
cial quantum systems over their solid state counterparts 
for studies of MBL. 

Artificial quantum systems containing relatively small 
number of quantum particles are also attractive from the¬ 
ory side as they can be simulated numerically by exact 
diagonalization of their Hamiltonians or by approximate 
methods. Recent theoretical work was focused on numer¬ 
ical and analytical studies of interacting one dimensional 
spins chains [8-18] that has shown that spin systems con¬ 
taining more than ten spins and involving thousands of 
many-body eigenstates exhibit the MBL behavior in suf¬ 
ficiently strong disorder. The MBL phase can be char¬ 
acterized by the existence of infinite number of local in¬ 
tegrals of motion [8-11], the entanglement [12-17], as 
well as the spectral properties of eigenstates [3]. While 
these characteristics can be measured in principle, the 
corresponding experiments are strenuous as they require 
a full quantum tomography of systems’ states. 

In this paper we propose an alternative strategy to 
identify the localization in a disordered system of inter¬ 
acting spins. We analyze the response of a pure state of 
a Heisenberg spin chain in random magnetic field along z 
axis to a sudden application of a magnetic field (quench) 
perpendicular to z that acts on a single spin. We evalu¬ 
ate the inverse participation ratio (IPR) of an eigenstate 
of unperturbed Hamiltonian in the basis of the quenched 
Hamiltonian. The IPR is small in the delocalized, or er- 
godic, regime when the initial state overlaps with many 
eigenstates of the new Hamiltonian. In the localized 
regime, an application of a quench does not affect ma¬ 
jority of eigenstates and the typical value of the IPR is 


about unity. Points at which the IPR starts increasing 
fast form a curve in the energy vs disorder strength plane 
and define the mobility edges. We note that due to finite 
size of our system, the mobility edge cannot be defined 
rigorously and the mobility edge has to be treated as a 
crossover region. Recently, the authors of Ref. [18] argued 
that even in infinite system MBL and ergodic phases are 
separated by a crossover region. 

Since the IPR is not easily measurable in experiments, 
we also investigate correlations in single spin measure¬ 
ments before and after the quench. The covariance be¬ 
tween these two measurements is small for delocalized 
states, but rapidly increases for localized states, as the 
quench only weakly affects configuration of far away 
spins. The mobility edge obtained from the covariance 
is consistent with the mobility edge obtained through 
the IPR, as well as through analysis of the entanglement 
entropy [19, 20]. Moreover, at strong disorder in the lo¬ 
calization regime, the spin response to the quench de¬ 
creases exponentially as a function of the distance of the 
monitored spin from the quenched spin. We utilize this 
exponential decay to evaluate the localization length as 
a function of disorder and demonstrate that the localiza¬ 
tion length exceeds the system size at the onset of MBL. 


i = L - 1 



FIG. 1. (Color online) Schematic picture of the one-di¬ 
mensional spin-1/2 chain with periodic boundary condition. 
Along the chain, each spin is subject to a random onsite field 
hi along 2 direction and the spins are coupled by nearest 
neighbor Heisenberg interactions with strength J. At time 
to when the local quench is turned on, a transverse magnetic 
field / is applied to one of the spin, labeled by i = 0. 

Local Quench.- We consider a Heisenberg chain of L 
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spins with random on-site field in the z direction and 
periodic boundary conditions: 

2—0 2—0 

where hi on each site is a random variable distributed 
uniformly in the interval [—Vb, W] and is the Pauli 
matrix for spin at site i, see Fig. 1. Throughout the 
paper, we use J as a fundamental unit and set J = 1. We 
denote eigenstates of Hq by |as^), where the total spin in 
2 direction Sz = '^ 2 *^ i® conserved for Hamiltonian (1). 

Previous numerical [17, 21] and analytical [22] studies 
were focused on the subspace defined by energy states 
|ao) with Sz = 0, where the MBL phase develops at 
W > 3.4. In this paper, we study the effect of the sudden 
quench V = fo-x~^\ see Fig. 1, applied to Hamiltonian 
(1) at site 1 = 0. The new Hamiltonian H = Hq + V 
break conservation of Sz, and system dynamics occur 
in the full 2^ dimensional space with the basis defined 
by eigenvectors \a) of H. The response of the system 
to this quench is reminiscent to the quasiparticle spec¬ 
tral function that appears in problems about transport 
in disordered electron systems with interactions. [23] 

We consider a system that was originally prepared as 
a pure state \'ijj) in the subspace of states with Sz |^) = 
0 and calculate its response over a long time after the 
quench: 

{O) = lim i r d{t)dt = Tr^ = Trpl^^j^O, (2) 

T-)-oo 1 Jq 

where p{t) = exp exp (iHt) is the density ma¬ 

trix and = p(t) is the time-averaged density matrix 
initialized in a pure state. Then, p(0) = ['0) (0|, 

Pde = K«IV')I^ Pa = \a){a\ (3) 

a. 

and Pa is the projection operator on new eigenstates jd) 
of H. The off-diagonal elements of vanish after time 

averaging and matrices belong to a diagonal ensem¬ 
ble [24]. 

Below we concentrate on initial states |0) that are 
eigenstates [ag) of the initial Hamiltonian Hq with Sz = 
0, here level index ag = 1,... N{L, 0) runs over eigenen- 
ergies of Hq ordered in increasing order, N{L,S) = 
(2^)!/[(2('^-‘5)/2)!(2('^+^)/2)!]. 

Inverse participation ratio - We demonstrate that the 
Hamiltonian H with the quench exhibits the MBL phase 
by analyzing the IPR: 

IPR„„ = ^ |(ag|d)|4 = Tr{pi3““>P„„} , (4) 

a=l 

where Pag = jao) {o^ol- The IPR is a measure of portion 
of the Hilbert space explored by the system after the per¬ 
turbation V is turned on and was introduced earlier for 


analysis of mobility edge in disordered electron systems 
with interactions [23, 25]. At weak disorder, the motion 
of the system is ergodic and the state spreads over a large 
fraction of the Hilbert space. As a result, the IPR is small 
~ 2“^. On the other hand, in the strong disorder limit, 
the ergodicity breaks down and the evolution of many- 
body wavefunctions is restricted to a small portion of the 
Hilbert space defined by the local integral of motion [8] . 

To investigate behavior of the IPR across the MBL 
transition we performed exact diagonalization for L = 12 
spins and N = 1000 realizations to obtain all eigenstates 
jag) for Hq and [d) for H. The quench strength is fixed 
at / = J throughout our computations. The average 
IPR as a function of disorder strength W and place of 
eigenenergy within the energy band e = Qfg/iV(L,0) is 
plotted in Fig. 2(a). After disorder averaging IPR(IF,e) 
clearly reveals the existence of a mobility edge that dis¬ 
tinguishes ergodic from localized states. To justify the 
nature of the mobility edge, we plot the histogram of the 
distribution of log 2 IPR in Fig. 2(b). In the weak and 
strong disorder limit, the distributions of log 2 lPR are 
highly concentrated at either oc {—L) or 0, respectively. 
However, log 2 IPR in the crossover region is broadly dis¬ 
tributed between (—L) and 0 with its standard deviation 
oc L so that the standard deviation is peaked at the mo¬ 
bility edge. This singular behavior is similar to that of 
the fluctuations of the entanglement entropy [17] and the 
form of the IPR distribution at weak and strong disorder 
resembles the distribution of the quasiparticle spectral 
functions [1, 23]. 

In the lower panels Fig. 2(c), we make two vertical 
cuts at fixed disorder strengths IT = 3, 7 on the phase 
diagram. For moderate disorder strength W = 3 and in 
the presence of mobility edge, the IPR reaches ~ 0.4 at 
the edges in an unsymmetrical manner, but it sharply 
drops to ~ 10“^ and forms a fiat basin of ergodic states 
in the middle of the band. At strong disorder IT = 7, 
IPRao ^ 0.4, indicating that most states are localized. 
When we fix the energy of an initial state in the middle 
of the band, e ~ 0.5, the standard deviation of IPR has 
a peak at the critical value ~ 3.5, when the upper 
and lower mobility edges merge, see Fig. 2(d). We note 
that at strong disorder, IT ~ 10, the standard deviation 
of IPR remain of the order of unity, see Fig. 2(b,d). in¬ 
dicating a broad distribution of the IPR in the localized 
regime. This behavior of the IPR in the middle of the 
energy band can be connected to the statistical orthogo¬ 
nality catastrophe [26, 27]. 

Local observables - While the above approach to de¬ 
tect mobility edge through IPR is suitable for numerical 
calculations, it is hardly realized experimentally. For ex¬ 
perimental detection of the mobility edge, we propose a 
different approach. We study the dynamics of single spins 
and investigate the correlation of their expectation values 
(ni*^) before and after a local perturbation. Our moti¬ 
vation is based on the previous observations [28] that in 
the MBL regime the eigenstate thermalization hypoth¬ 
esis (ETH) is violated and information about the local 
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FIG. 3. (Color online) Scatter plot of (Pfa), QJa)) (a) 
the localized regime with VF = 10 and (b) the ergodic regime 
with W = 2.5 for states |q;o) in the middle of the band. Data 
obtained with L = 12 and N = 1000. The corresponding 
probability distributions of Pi and Qi for the localized regime 
(c) in which the single spin measurement with and without 
the quench are similar and for the ergodic regime (d) the 
discrepancy in the distribution with and without the quench 
indicates thermalization. 


FIG. 2. (Golor online) (a) Average IPR over 1000 disorder 
realizations as functions of disorder strength W and energy 
density e for a system of L = 12 and / = J = 1. The 
many body mobility edge along the outer border of the light- 
gray sector encloses a region of ergodic states with IPR 
2“^. The horizontal and vertical lines indicate parameters 
for panels (c) and (d). (b) The histogram of logj IPR in the 
middle of the band for W = 0.5, 4 and 10. The distributions 
are narrow for weak and strong disorder, but broad near the 
critical region, (c) The IPR at fixed disorder W = 3, when the 
states in the middle of the band are ergodic but localized near 
the band edges, and W — 7 when all states are localized, (d) 
The IPR and its standard deviation for fixed energy density 
e~0.5. 

integrals of motion maintains correlations between spin 
states before and after the quench. Otherwise, at weak 
disorder the motion is ergodic and information about ini¬ 
tial states is lost. 

For an interacting system, the observables set by a fi¬ 
nite degrees of freedom can be evaluated by the reduced 
density matrices in which the off-diagonal elements are 
essentially zero due to dephasing even if the system starts 
in some arbitrary pure state. Alternatively, time aver¬ 
aged expectation values of a local operator are also char¬ 
acterized by a diagonal ensemble of the density matrix 
see Eq. (2). One can measure such local observables 
before and after the quench. Below, we again consider 
a special case of initial states that are eigenstates |ao) 
of Hq with = 0, although the results remain qualita¬ 
tively the same for a system that is initially prepared in 
an arbitrary pure state. We calculate expectation values 


of of spin i before and after the quench at site 0: 

Ko =io‘o\cri"Uao): (5a) 

Qao l(«o|d)|^ (a|crW ^ (5b) 

a. 

where Pa |(d|ao)|^ 

A good distinction between ergodic and MBL phases 
can be obtained by analyzing correlations between mea¬ 
surements of and Qag- Running over all eigenstates 

|ao) in the = 0 sector, we collect P^^ and for 
a number of disorder realizations and present the scat¬ 
ter plot for pairs of {P^^, Qh^) of spin z = 0 (directly 
perturbed spin, open circles) and i = L/2 — 1 (the far¬ 
thest spin from the quench, filled squares) in the mid¬ 
dle of the band, ao ~ N{L, 0)/2 for strong {W = 10) 
and weak {W = 2.5) disorder, see Fig. 3. In the local¬ 
ized phase at strong disorder, the eigenstates are prod¬ 
uct states consisting of physical spins |q;) = |4, (t))^ 

and therefore the local spin projection is good quantum 
number P^,^ ~ ±1. Provided that the quench is smaller 
than the local random field, the eigenstates of i?o and 
H differ by terms 0{f jW), resulting in an almost un¬ 
changed ~ when the distributions of P^^ and 
are peaked at ±1, see Fig. 3(c). On the scatter 
plot, pairs (Pi'/^“\ are distributed along the 

line Pal"^ ^ = Qao ^ ^, indicating the two measurements 
are strongly correlated. In the ergodic phase at weak dis¬ 
order {W = 2.5), the eigenstate is a linear combination 
of a large number of product states, and therefore local 
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spin projection is not a good quantum number. Two sets 
of local spin measurements form an elliptic cloud, see 
Fig. 3(b). This distribution points to ergodicity: the dis¬ 
tribution of single spin measurement for all possible spin 
configurations have maxima at = 0 and = 0 
and the distribution of QLq narrower due to the equi¬ 
libration of the system after the quench, see Fig. 3(d) for 
W = 2.5. Remarkably, the distributions are almost indis¬ 
tinguishable for MBL regime at W = 10, cf. Figs. 3(c,d). 

Correlations between P^^ and are characterized 
by the covariance with respect to disorder realizations: 

Cap (i) = iP^„ -Pao)' (Qao -Qao) = 

P^^ ~ Qap — 0- We use the covariance to map out 
the phase diagram as a function W and e = ao/iV(L,0), 
as shown in Fig. 4. In the ergodic regime the averaged 
value C'ao(T/2) vanishes but it saturates to 1 deep in the 
MBL phase where both P^^ and take almost identi¬ 
cal values (±1). Similar to the IPR, Cag{i) reveals the 
many-body mobility edge, marked by the border of the 
black region in Fig. 4(a). 



FIG. 4. (Color online) (a) Correlator ^ 

function of disorder W and energy density e for L = 12 and 
2000 realization. The many-body mobility edge is marked in 
white, (b) D' , see eq. (6) as a function of distance between the 
monitored and quenched spins (points), and the correspond¬ 
ing fitting curve, Eq. (7) (dashed lines), (c) The many-body 
localization length ^ extracted from Dag as a function of dis¬ 
order strength in the middle of the band, qq — N{L,Q)/2. 
The dashed line corresponds to ^ = 1 to show that the at 
strong disorder the system is localized on the atomic scale. 


With decreasing the spatial separation between the 
quenched spin and the monitored spin, we observe larger 
deviations between P^^ and QLp perturbed spin 

can be thermalized with subsystem of surrounding spins 
within a cluster of a size of the localization length. On 
the contrary, in the ergodic regime, the distribution of 
{Pag, Qag) is insensitive to the spatial separations be¬ 


tween quench at site 0 and monitored spin f, suggesting 
that the thermolization happens throughout the whole 
system. The ensemble statistics of and Qag can be 
used to evaluate the localization length. Intuitively the 
localization length is a scale at which spin texture form 
localized clusters and thereby the ergodicity is broken. 
Essentially deep in the localized regime, localization oc¬ 
curs on atomic scales with ^ 1. With decreasing dis¬ 

order, the localization length grows and once the scale 
is beyond the system size, the entire system cannot be 
decomposed into independent clusters and the ergodic 
behavior takes place. Therefore, the localization length 
is an indicator of the onset of MBL regime that can be 
determined by the spatial sensitivity of the response to 
the local quench. The deviation between the measure¬ 
ment Pa^ and Q’ag averaged over disorder realizations is 
given by the “distance”: 


D\ = ^{Pi^g-QlgY- ( 6 ) 

We argue that in the localized regime Da^ is an expo¬ 
nentially decaying function with respect to the distance 
between spin i to the quenched spin: 

D\glD\ cs e-*/« + (7) 

where ^ is the localization length. The first two expo¬ 
nential terms in Eq. (7) arise from the periodic boundary 
condition, and the last term is a normalization constant. 
In Fig. 4(b) we present the ratio D'’ jD'^ and its best fit 
with respect to choice of ^ in Eq. (7) as a function of spa¬ 
tial separation i for several different disorder strengths 
W = 3, 5,10 in the middle of the band. In Fig. 4(c), we 
illustrate the fitted localization length ^ for different val¬ 
ues of disorder W. Due to the finite size of the system, 
at the critical disorder Wg ~ 3.4 the localization length 
^ does not diverge but becomes larger than the system 
size, consistent with our argument about the crossover 
between ergodic and MBL phases. On the other hand, 
^ saturates to unity at strong disorder deep in the MBL 
phase. 

Summary - We showed that as a result of local quench, 
the absence of thermalization in the MBL phase can be 
characterized by the IPR of the eigenstates of the orig¬ 
inal Hamiltonian in the basis of the quenched Hamilto¬ 
nian. In particular, the IPR fluctuations are enhanced 
at the crossover between the ergodic and MBL phases. 
Our analysis of the single spin measurements in response 
to a quench shows a plausible experimental technique to 
search for the mobility edge and the localization length. 
This approach does not require substantial measurements 
of a quantum system as compared to more complicated 
measurement of many-particle entanglement. Our anal¬ 
ysis demonstrated that for a system with MBL behav¬ 
ior, simple single spin measurement can reveal the in¬ 
dispensable characteristics, complementary to more so¬ 
phisticated routes based on the growth of entanglement 
entropy [12, 13] or quantum revivals[29]. 
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